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We consider the Gross-Neveu model with a continuous chiral symmetry in two and three space- 
time dimensions at zero and finite temperature. In order to study long-range phase coherence, we 
construct an effective low-energy Lagrangian for the phase 9. This effective Lagrangian is used to 
show that the fermionic two-particle correlation function at finite temperature decays algebraically 
in 2+1 dimensions and exponentially in 1+1 dimensions. 
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Introduction. Low-dimensional field theories have 
been studied extensively as toy models for QCD at zero 
and finite temperature. Many of these theories are re- 
markably rich, sharing a number of properties with four- 
dimensional QCD. For example, both the Gross-Neveu 
model and the O(N) nonlinear sigma model in 1+1 di- 
mensions are asymptotically free and have nonperturba- 
tively generated mass gaps. In addition, the latter has 
instanton solutions for N = 3. In contrast with QCD in 
3+1 dimensions, these models are renormalizable in the 
X/N expansion both in 1+1 and 2+1 dimensions. Thus 
one can systematically study these models in a nonper- 
turbative setting 0, Q. There is, however, an impor- 
tant difference between QCD and low-dimensional field 
theories, namely the possibility of spontaneous break- 
ing of continuous global symmetries. In QCD with Nf 
massless fermions, chiral symmetry is spontaneously bro- 
ken at zero temperature and this gives rise to Nj — 1 
massless Goldstone particles. For Nt = 2, these are 
the well known pions. Chiral symmetry is expected to 
be restored at a critical temperature of approximately 
150 MeV depending on the number of flavors. In 1+1 
and 2+1 dimensions, the Mermin- Wagner theorem 0,Q 
forbids the spontaneous breakdown of continuous sym- 
metries at any finite temperature . The reason is 
that the phase fluctuations are so strong that they de- 
stroy the presence of a condensate. The role of phase 
fluctuations was studied in the seminal papers on the 
Kosterlitz-Thouless phase transition of vortex unbind- 
ing 0, 0] . The low-temperature phase where the vortices 
are bound in pairs is characterized by an algebraic de- 
cay of the one-particle density matrix and superfluidity. 
The high-temperature phase of unbound vortices is char- 
acterized by loss of superfluidity and exponential fall-off 
of correlations. More recently, phase fluctuations have 
been studied in connection with Bose-Einstein condensa- 
tion of atomic gases in harmonic traps in one -and two- 
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1 In 2+1 dimensions, spontaneous breaking of continuous symme- 
tries may occur at T = 0. Consequently, the critial temperature 
for the restoration of a continuous symmetry is T c = 0. 



dimensional Bose gases 0, EJ EH- In the following, 
we consider the Gross-Neveu model with a continuous 
chiral symmetry |llj. Many of its properties have been 
examined in detaiT|l HJ El 13 El El In the present 
paper we discuss the possible phases of the theory. By 
using density and phase variables, one can conveniently 
address the issues of long-range order and phase struc- 
ture as a function of temperature. These issues were 
discussed in Refs. ^2 E^ : nere we present a somewhat 
more general treatment at zero and finite temperature. 

Gross-Neveu model. The Euclidean Lagran gian of 
the Gross-Neveu model with a U(l) symmetry is [ll| 
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where j — 1,2, „JV. In 2+1 dimensions, the 7-matrices 
are 4x4 matrices and in 1+1 dimensions they are 2 x 
2 matrices Q- They satisfy {7^,7^} = 2<Jy and 75 = 
~ 27071 ■ The Lagrangian has a U(l) chiral symmetry: 

^ -> e^Vi , (2) 

where <fi is a constant phase. Following |ll|. we introduce 

2 - 2 - 

the auxiliary fields a = jfipj^j an d t = ^f'4>jil5' l l , j- The 
Lagrangian can then be written as 



tpj (cr + i7T7 5 ) tpj . (3) 



In terms of the fields a and 7T, chiral symmetry can be 
written as 



(cr + in) — ► (cr + in) e 2 ^ . 
Defining the density p and phase 6 variables by 

a + in = pe w , 
and the new fermion field \ by 



XL = e , XR 



- p if>/2 



iJ'I 



we can write the Lagrangian Q as 



N 



(4) 
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In terms of p and 9, the chiral symmetry is simply 9 —> 
9 + c, where c is a constant. In these variables the chiral 
symmetry is manifest. Moreover, the operator $ acting 
on 9 in (JJ ensures that the field is massless and that it is 
derivatively coupled. The auxiliary field p is now written 
as a sum of a space-time independent background po 2 
and a quantum fluctuating field p; p = po + p- Note 
that a nonzero expectation value for p does not violate 
chiral symmetry since p is invariant under global U(l) 
transformations. Since the Lagrangian 10 is quadratic 
in the fermion fields, we can integrate over \, and we 
obtain the effective action for p and 9: 



S = -ATrln 



N 



<j) + ^75^0 + po+p 
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where (3 = 1/T is the temperature (in units where h = 1) 
and d — 1 is the dimension of space. The trace is over 
both Dirac indices and spacetime. The next step is to 
expand the functional determinant around the classical 
solution of the equation of motion 9 — and p = 0. This 
expansion generates propagators for p and 9 as well as 
interaction vertices among them. 

At leading order in 1/N, the auxiliary fields p and 
9 do not propagate. However, at next-to-leading or- 
der, quantum fluctuations induce the inverse propaga- 
tors. These are obtained from expanding the functional 
determinant in JSJ to second order in the fluctuations: 



n P (Q) 
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n e (Q) = lTv(I)Np 2 jf 



(9) 
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(10) 

where Tr(J) is the trace of the identity matrix /. We 
have introduced the notation 
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P = (2n+1) 



jd-1 



P 
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(11) 



and P = {Po , p) is the Euclidean momentum. The inte- 
gral over spatial momenta is regularized by using an ul- 
traviolet cutoff A. Eq. © is ultraviolet divergent and the 
divergence can be eliminated by renormalizing the cou- 
pling constant g. By analytic continuation to Minkowski 



2 Note that in 2+1 dimensions, po is nonzero at T = only if 
the coupling g is larger than a critical coupling g c Q In 1+1 
dimensions, po is always nonzero at T = 0. 



space (Po — > iuj + iff) , one can examine the analytic 
structure of the propagators for p and 9. For example, 
at zero temperatu re, both p ropagators have a branch cut 
starting at oj = ^/q 2 + ^p\. The 9 propagator has in ad- 
dition a pole at w = q, which ensures that the field is 
massless. This pole remains at finite temperature, which 
follows immediately from Eq. (|10(l . 



The effective potential V through leading order in 
1/JVis 



V = 



Np 
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The value of po is determined by minizing the effective 
potential and to leading order in the 1 /N expansion this 
leads to the gap equation 



1 



Tr(/)rf 
7i 



1 



\p } P 2 + Pt 



(13) 



This gap equation coincides with the gap equation 
in Gross-Neveu model with a discrete Z2 symmetry 
Eq. (|13|l is ultraviolet divergent and it is rendered finite 
by renormalizing the coupling constant g. At T = 0, the 
gap equation can be solved explicitly for po- For exam- 
ple, in 1+1 dimension, we obtain p\ = A 2 e~ 27r / 9 «, where 
gu is the renormalized coupling and A is the ultraviolet 
cutoff. Finally, note that one can also use the gap equa- 
tion i|13[l to eliminate the divergent sum-integral in the 
expression for the p-propagator. 

If we are interested in the long-distance properties 
of the theory Q), we can simplify calculations by notic- 
ing that for momenta p much smaller than po, the mas- 
sive field p decouples and we are left with an effective 
long-distance theory for the phase 9. To leading order in 
derivatives, this can be written as a free field theory: 



\*(d»9f 



(14) 



where a is the stiffness of the phase fluctuations and is 
given by 



Tr(/)Ap 2 ^ 



1 



(P 2 + P 2 ) 
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(15) 



The summation in (|15fl is over the bosonic Matsubara 
modes Po = ^nnT. The expression for a follows imme- 
diately from i|l(J[l . Chiral symmetry, 9 — > 9 + c, forbids 
terms in the effective Lagrangian not involving deriva- 
tives of 9 such as 9 4 and 9 2 (V9) 2 . Thus the next opera- 
tor in the low-energy expansion involves four derivatives 
and the results concerning the infrared properties of the 
Gross-Neveu model obtained using (|14|) are stable against 
perturbations [l2|. 

In order to discuss the possible phases of the Gross- 
Neveu model, we must examine the long-distance behav- 
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ior of the following four-fermion correlation function [T^| 

p 2 (x,0) = (^(l + TsW^WCl-TlsWO)) ■ 



(16) 



The function p2(x,0) can be written as 

p 2 (x,0) = (p(x)e-^p(0y e W 

Ho)] 2 ) 



(17) 



where we in the second line have replaced p by its ex- 
pectation value po and used Wick's theorem to rewrite 
the exponent. This replacement is justified since den- 
sity fluctuations are negligible for long distances. Using 
the effective Lagrangian (|14|) . the exponent (|18fl can be 
written as 



({8(x)-e(o)f) 



1 r dp- 1 

a J (^iy- 1 
fl — cos(p • x)] 



P t-u ( @P 
coth 



x 



P 



(18) 



The behavior of (|18fl in the limit x — > oo then determines 
the possible phases of the Gross-Neveu model. 

2+1 dimensions. We first notice that the gap equa- 
tion (|13fl and the stiffness (|15fl can be calculated explic- 
itly: 



-1 = -£»-?Iln[l 
9r it it l 

a = — p tanh— . 



-l3po] 



(19) 
(20) 



Averaging over the angle between p and x in (|18fl . we 
obtain 



(ie(x)-e(o)f) 



i 



dp coth 



0p 



1 



J (px) 
(21) 



where Jo(px) is a Bessel function of the first kind. At 
T = 0, the expression (|2~TJl reduces to 



(i9(x)-9(0)f) = - 



2irax 



(22) 



In the large- a: limit, the exponent vanishes and this shows 
that there is a real condensate for T = 0. This is the 
usual phase with the U(l) symmetry spontaneously bro- 
ken and long-range order 0. At finite temperature the 
exponent l|21|l can not be evaluated in closed form. How- 
ever, one can expand the integrand and integrate term 
by term. This gives the following convergent series 



([0(x)-9(O)f) 
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The large-x behavior can be found by integrating the sec- 
ond term in (|23|) considering n as a continuous variable. 
Alternatively, one can the apply dimensional-reduction 
argument that for distance scales x much larger than 
1/T and l/po, the correlator is dominated by the zeroth 
Matsubara mode. Either way, in the limit x — ► oo one 
obtains 



([9(x)~9(0)f) = L. 

(XK 



1 X 



C 



(24) 



where C is a constant. The correlator decreases alge- 
braically with exponent 7 = T/2ira. This shows that 
at finite temperature there is a quasicondensate whose 
density is pq. In this phase, there is a local gap, but the 
phase is incoherent. This is referred to as a pseudogap. 

The gap equation 1)13(1 ceases to have a nontrivial 
solution for On for temperatures T above T p = po(T = 
0)/21n2 dEi ( Note tnat t p is the critical temperature 
for the phase transition in the case of the 2+1 dimen- 
sional Gross-Neveu model with a Z 2 -symmetry) . Above 
this temperature, the leading-order effective potential 
then reduces to the ideal-gas value, V = — 3T 3 £(3)/47r. 

1+1 dimensions. Averaging over the angle between 
p and x in l|18(l . one obtains 

<[*(*) -*(0)]») = ^fVoth(^ 
[1 — cos(px)] 



p 



(25) 



At zero temperature and in the large- a: limit, the expo- 
nent reduces to 



({9(x)-9(0)f) = —Inx 



(26) 



At zero temperature, the stiffness (|15|) reduces to a = 
N/2n which implies that p 2 {x, 0) = p\x~ x l N . This agrees 
with Refs. Thus at T = 0, the two-particle cor- 

relation function decays algebraically and one has "al- 
most long-range order" . At finite temperature the expo- 
nent I125H can be calculated exactly and reads 



([9(x)-9(0)f) 



1 



2na 

where K is a constant. We then obtain 

1 

2a 



[lnsinh(7rav73) + K] , (27) 



lim {[9{x) - 0(O)r) 



-Tx 



(28) 



which shows that the correlator decays exponentially, 



p 2 (z,0) 



e 



-Tx 



Thus not even a quasicondensate is 



present at finite temperature. The large-x (|28|l behavior 
could also have been derived by calculating the contribu- 
tion from the static Matsubara mode, which dominates 
in this limit. 

For temperatures above T p — po(T — 0)e 7 - B /tt [l5| . 
the gap equation has only the trivial solution po — (The 



3 



critical temperature for the phase transition in the dis- 
crete case is T c = 0). The leading-order effective poten- 
tial therefore reduces to the ideal-gas value, V = — nT 2 /6. 
More generally, due to asymptotic freedom we expect in 
the high-temperature limit that the 1/N expansion re- 
duces to perturbation theory. 

Summary. In the present paper, we have studied the 
phase fluctuations of the [/(l)-symmetric Gross-Neveu 
model in 1+1 and 2+1 dimensions. Constructing an ef- 



fective low- momentum field theory for the phase 9, we 
have calculated the two-particle correlation function. In 
2+1 dimensions, there is a real gap at zero temperature, 
while the phase fluctuations turn it into a pseudo-gap 
at finite temperature. In 1+1 dimensions, the quasicon- 
densate which is present at T = 0, is destroyed at finite 
temperature due to phase fluctuations. In particular, this 
implies that the Kosterlitz-Thouless-like phase present at 
T = |12| disappears at finite temperature. 
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